The Bose-Einstein correlation function C^Q) from a Quantum Field Theory point of 

view 
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We show that a recently proposed derivation of Bose-Einstein correlations (BEC) by means of a 
specific version of thermal Quantum Field Theory (QFT), supplemented by operator-field evolution 
of the Langevin type, allows for a deeper understanding of the possible coherent behaviour of the 
emitting source and a clear identification of the origin of the observed shape of the BEC function 
C*2((5). Previous conjectures in this matter obtained by other approaches are confirmed and have 
received complementary explanation. 
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In this work we would like to focus attention on 
two specific features of Bose-Einstein correlations (BEC) 
clearly visible when BEC are presented in the language 
of some specific (thermal) version of Quantum Field The- 
ory (QFT) supplemented by the operator-field evolution 
of Langevin type proposed recently [1,2]. Because the 
importance of BEC and their present experimental and 
theoretical status are widely known and well documented 
(see, for example, [3] and references therein), we shall 
not repeat it here. The features we shall discuss are: 
(i) how the possible coherence of the hadronizing system 
(modelled here by some external stationary force occur- 
ring in the Langevin equations describing hadronization 
process [2]) influences the BEC function C*2(Q) and (n) 
what is the true origin of the experimentally observed 
Q-dependence of the C2(Q) correlation function in the 
approach used here . 

In what concerns the first point we have obtained iden- 
tical expression as derived in [4] by means of a quantum 
statistical (QS) approach with novel interpretation of the 
chaoticity parameter p introduced there. On the other 
hand, our results (and that of [4]) differ from the formula 
recently obtained in [5]. We shall argue that the origin 
of this difference lies in different ways of introducing the 
concept of coherence in both approaches. 

In the second point we demonstrate that in order to 
obtain a given (experimentally observed) shape of the 
BEC correlation function C2(Q) (i.e., its Q-dependence, 

where Q — |fc M — k'^\ — (fc M — fc^) 2 ) one has to ac- 
count somehow for the finiteness of the space-time re- 
gion of the particle production (i.e., of the hadronizing 
source). In our approach this means the necessity of 
smearing out of some generalized functions (delta func- 



tions: <5(<3 A 



k')) appearing in the definition of 



thermal averages of some operators, which is character- 
istic feature of QFT approach used here. The freedom 
in using different types of smearing functions to perform 
such a procedure allows us to account for all possible dif- 
ferent shapes of hadronizing sources apparently observed 
by experiment. [6]. 

Referring to [1] for details, let us recapitulate here the 



main points of our approach. The collision process pro- 
duces a lot of particles out of which we select one (we as- 
sume for simplicity that we are dealing only with identical 
bosons) and describe it by operator b(k, t) (the notation 
is the usual one: b(k, t) is an annihilation operator, k is 
3- momentum and t is a real time). The rest of the parti- 
cles are then assumed to form a kind of heat bath, which 
remains in equilibrium characterized by a temperature 
T — 1//3 (which will be one of our parameters). All av- 
erages ((. . .)) are therefore thermal averages of the type: 
((. . .)) = Tr [(. . .)e-P H ] /Tr (e"^). We shall also al- 
low for some external (to the above heat bath) influence 
to our system. Therefore we shall represent the operator 
b(k, t) as consisting of a part corresponding to the action 
of the heat bath, a(k,t), and also of a part describing 
action of these external factors, R(k,t): 



b(k,t) = a{k,t) +R(k,t). 



(1) 



The time evolution of such a system is then assumed to 
be given by a Langevin equation [2] 



id t b(k, t) = F{k, t) - A{k, t) + P 



(2) 



(and a similar conjugate equation for b + (k,t)). These 
equations are supposed to model all aspects of the 
hadronization process. The combination F{k, t) — A(k, t) 
represents the so called Langevin force and is therefore 
responsible for the internal dynamics of hadronization in 
the following manner: A is related to stochastic dissipa- 
tive forces and is given by [2,1] 



A{k,t) 



4-oo 



drK(k,t-T)b(k,T), 



(3) 



with the operator K {k, t) being a random evolution field 
operator describing the random noise and satisfying the 
usual correlation-fluctuation relation for the Gaussian 
noise [8]. The operator F(k,t) describes the influence 
of heat bath, 



F(k,t) 



'dw 
2^ 



tp(k M )c(k l j)e 



(4) 
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Our heat bath is represented by an ensemble of damped 
oscillators, each described by operator c(fc M ) such that 
[c(k tl ),c + (k' ll )] = S 4 (kf l — k' ), and characterized by some 
function V(fc At ) [9]. Finally, the constant term P (repre- 
senting external source term in Langcvin equation) de- 
notes the possible influence of some external force (as- 
sumed here to be constant in time). This force would re- 
sult, for example, in a strong ordering of phases, leading 
therefore to the coherence effect in the sense discussed in 
[4]. Out of many details (for which we refer to [1]) what 
is important in our case is the fact that the 2-particle 
correlation function for like-charge particles, is defined 
as (hp = (tjj = k°, kj)): 



C 2 (Q) 



f{kfj,, k') 



/(V 

[1 + 



D(k„k'^} 



(5) 



where RfyK) _ = (b+(k^b+(k'^b(k^b(k^) and 
f(kfi) = {b + {k fl )b{k fl )) are the corresponding thermal 
statistical averages (in which temperature T enters as a 
parameter) with b(k^) = a(k^) + R(kn) being the corre- 
sponding Fourier transformed stationary solution of eq. 
(2). As shown in [1] (notice that operators R(kp) by 
definition commute with themselves and with any other 
operator considered here): 



+ <a+(fe M )o(^))<a + (^)a(fe M )) + 

+ (a+(k^a(k'^)R+(k^R(k^ + 
+ (a + (k^a(k^)R + (k^R(k^), 



/(*;„) = (~a+(k,)a(k,)) + \R(k„)\ 2 . 



(6) 
(7) 



This defines D(fc„, = /(*;„, ■ Rk'J] - 1 in 

(5) in terms of the operators a(fc M ) and i?(fc Al ), which in 
our case are equal to: 



a(fe M ) 



F{kf * ) - and R{k„) 



K{k») 



K{k») 



- (8) 



The multiplicity N depending factor £ is in our case equal 
to £(N) = (N)' 2 /{N(N - 1)). This means therefore that 
the correlation function C%{Q), as defined by eq. (5), 
is essentially given in terms of P and the two following 
thermal averages for the F(k,t) operators: 



(F+(k,t)F(k',t')) =6 3 (k-k') 
(F(k,t)F+(k',t')) =S 3 (k-k' 



g n( w )e*-<*-0 (9) 



/ 



du> 
2^ 



|^| 2 [1 + n(w)]e- iw (*-*') 



where n(ui) — {exp [(ui — fi)(3] — 1} 1 is the number of 
(by assumption - only bosonic in our case) damped os- 
cillators of energy oj in our reservoir characterized by 
parameters \i (chemical potential) and inverse tempera- 
ture (i = 1/T (both being free parameters) [10]. Notice 
that with only delta functions present in (9) one would 
have a situation in which our hadronizing system would 
be described by some kind of white noise only. The in- 
tegrals multiplying these delta functions and depending 
on (a) momentum characteristic of our heat bath ^(k^) 
(representing in our case, by definition, the hadronizing 
system) and (b) assumed bosonic statistics of produced 
secondaries resulting in factors n{u>) and 1 + n(u>), re- 
spectively, bring the description of our system closer to 
reality. 

It should be stressed at this point that, contrary to 
the majority of discussions of BEC [3,4], we are working 
here directly in phase space [11], so far no space-time 
considerations were used. It is easy to realize now that 
the existence of BEC, i.e., the fact that C^Q) > 1, is 
strictly connected with nonzero values of the thermal av- 
erages (9). However, in the form presented there, they 
differ from zero only at one point, namely for Q — (i.e., 
for fcp = k'p). Actually, this is the price one pays for the 
QFT assumptions tacitly made here, namely for the infi- 
nite spatial extension and for the uniformity of our reser- 
voir. But we know from the experiment [3] that C2(Q) 
reaches its maximum at Q = and falls down towards its 
asymptotic value of Ci = 1 at large of Q (actually already 
at Q ~ 1 GeV/c). To reproduce the same behaviour by 
means of our approach here, one has to replace delta 
functions in eq. (9) by functions with supports larger 
than limited to a one point only. This means that such 
functions should not be infinite at = k^ — k'^ = but 
remain more or less sharply peaked at this point, oth- 
erwise remaining finite and falling to zero at small, but 
finite, values of |<5 M | (actually the same as those at which 
C2(Q) reaches unity): 



8{k^ - fcp =► n ■ y/Q(q = Q. r ). 



(11) 



(10) 



Here fio has the same dimension as the 5 function (ac- 
tually, it is nothing else but 4-dimcnsional volume re- 
stricting the space-time region of particle production) 
and 0(g) is a dimensionless smearing function which con- 
tains the g-dependence we shall be interested in here. 
In this way we are tacitly introducing a new parameter, 
r M , a 4- vector such that (r M ) 2 has dimension of length 
and which makes the product Q ■ r = Q^r^ = q dimen- 
sionless. This defines the region of nonvanishing density 
of oscillators c, which we shall identify with the space- 
time extensions of the hadronizing source. The expres- 
sion (11) has to be understood in a symbolic sense, i.e., 
that Q(Q ■ r) is a function which in the limit of r — > oo 
becomes strictly a S function. Making such replacement 
in eq. (9) one must also decide how to accordingly adjust 
n(u>) occurring there because now, in general, lu ^ u>'. In 
what follows we shall simply replace n(w) — > n(w) with 



2 



uj = (w + lj')/2 (which, for classical particles would mean 
that n(uj) — > y/ n(uj)n{u)')) . 

In such way r becomes new (and from the phenomeno- 
logical point of view also the most important) parameter 
entering here together with the whole function fl(Q ■ r), 
to be deduced from comparison with experimental data 
[12]. With such a replacement one now has 



(1 + <*)(! + a') 



(12) 



where 
Cl(q)= 7 -Q(q), 7 



n 2 (u>) 



a oc 



(13) 



with n(u>) the same as defined above. The parameter 
a is another very important parameter, which summa- 
rizes our knowledge of other than space-time characteris- 
tics of the hadronizing source (given by f2(g) introduced 
above). In particular it contains the external static force 
P present in the evolution equation (2). It is combined 
(in multiplicative way) with information on the momen- 
tum dependence of the reservoir (via ^(k^)] 2 ) and on 
the single particle distributions of the produced particles 
(via n(u> = fixcoahy) where [It and y are, respectively, 
the transverse mass and rapidity). Notice that a > 
only when P^O. Actually, for a = one has 



1 < C 2 (Q) < l+7ft(Q-r), 



(14) 



i.e., it is contained between limits corresponding to very 
large (lower limit) and very small (upper limit) values of 
P. Because of this a plays the role of the coherence pa- 
rameter [3,4]. For 7 ~ 1, neglecting the possible energy- 
momentum dependence of a and assuming that a' = a 
one gets the expression 



C 2 (Q) 



2a 



(1 + a) 2 



■V%) 



1 



(1 + a) 2 



Sl(q), (15) 



which is formally identical with what has been obtained 
in [4] by means of QS approach. It has precisely the same 
form, consisting two Q— dependent terms containing the 
information on the shape of the source, one being the 
square of the other, each multiplied by some combination 
of the chaoticity parameter p = 1/(1 + a) (however, in [4] 
p is defined as the ratio of the mean multiplicity of parti- 
cles produced by the so called chaotic component of the 
source to the mean total multiplicity, p = (N c h)/(N)). 
In fact, because in general a ^ a' (due to the fact that 
uo ^ uj' and therefore the number of states, identified here 
with the number of particles with given energy, n(uj), are 
also different) one should rather use the general form (5) 
for C 2 with details given by (12) and (13) and with a de- 
pending on such characteristics of the production process 
as temperature T and chemical potential p occurring in 
definition of n(w). 



Notice that eq. (15) differs from the usual empirical 
parameterization of C 2 (Q) [3], 



C 2 (Q) = l + A- fi(Q-r), 



(16) 



with < A < 1 being a free parameter adjusting the 
observed value of C 2 (Q = 0), which is customary called 
"incoherence", and with il(Q ■ r) represented usually as 
Gaussian. Recently eq. (16) has found strong theoretical 
support expressed in great detail in [5] . The natural ques- 
tion arises: which of the two formulas presented here is 
correct? The answer is: both are right in their own way. 
This is because each of them is based on different ways 
of defining coherence of the source. In [5] one uses the 
notion of coherently and chaotically produced particles 
or, in other words, one divides hadronizing source into 
coherent and chaotic subsources. In [4] one introduces in- 
stead the notion of partially coherent fields representing 
produced particles, i.e., one has only one source, which 
produces partially coherent fields. Our approach is simi- 
lar as we describe our particle by operator b(k, t), which 
consists of two parts, cf. eq. (1), one of which depends 
on the external static force P. The action of this force is 
to order phases of particles in our source (represented by 
the heat bath). The strength of this ordering depends on 
the value of the external force P. In any case, for P ^ 0, 
it demonstrates itself as a partial coherence [13]. 

Let us return to the problem of Q-dependence of BEC. 
One more remark is in order here. The problem with 
the 8{k tl — k! ) function encountered in two particle dis- 
tributions does not exist in the single particle distribu- 
tions, which arc in our case given by eq. (7) and which 
can be written as f(k fl ) oc (a + (fc M )a(fc M )) + ^(fc^)! 2 ~ 
(1 + a)(a + (k t _ l )a(k t _ l )) [14]. To be more precise 



/(*;„) = (l + a)-S(fc M ,Mi 



(17) 



where ^(/c^, fc M ) is one-particle distribution function for 
the "free" (undistortcd) operator a(k^) equal to 



kft) = n Q 



n(oj). 



(18) 



Notice that the actual shape of /(fc M ) is dictated both 
by n(uo) = n(u;T,fi) (calculated for fixed temperature 
T and chemical potential /i at energy uj as given by the 
Fourier transform of random field operator K and by 
shape of the reservoir in the momentum space provided 
by ip(k^)) and by external force P in parameter a. They 
are both unknown, but because these details do not en- 
ter the BEC function C 2 (Q), we shall not pursue this 
problem further. What is important for us at the mo- 
ment is that both the coherent and the incoherent part 
of the source have the same energy-momentum depen- 
dence (whereas in other approaches mentioned here they 
were usually assumed to be different). On the other hand 
it is clear from (17) that (N) = (N ch ) + (N coh ) (where 
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(N c h) and (N co h) denote multiplicities of particles pro- 
duced chaotically and coherently, respectively) therefore 
justifying definition of chaoticity p mentioned above. 

For an illustration we plotted in Fig. 1 the corre- 
lation function C^Q) as given by eq. (15) for differ- 
ent choices of £l(q) corresponding to different hadroniz- 
ing sources discussed in the literature [15-17] (here r = 

• Gaussian: Cl(q) = cxp (— Q 2 r 2 ) ; 

• exponential: 0,(q) = exp(— Qr); 

• Lorentzian: £l(q) = 1/ (1 + Qr) 2 ; 

• given by Bessel function [16]: 

All curves are drawn for the same values of the size pa- 
rameter r = 1 frn and assuming for simplicity constant 
and equal values of a and a' parameters (i.e., using sim- 
plified eq. (15)), which have been put equal a = 0.2 
here, just for illustrational purpose (it corresponds to 
p = 0.8 in [4]). Fig. 2 shows in detail (using Gaus- 
sian shape of fi(<z) function) the dependence of C%(Q) on 
different values of a = 0, 0.25, 1, 4 (again, used in the 
same approximate way as before and corresponding to 
p = 1., 0.8, 0.5, 0.2) and compare it to the case when 
the second term in eq. (15) is neglected, as is the case in 
majority of phenomenological fits to data. 

To summarize: using a specific version of QFT supple- 
mented by Langevin evolution equation (2) to describe 
hadronization process [1,2] we have derived the usual 
BEC correlation function in the form explicitly showing 
the origin of both the so called coherence (and how it 
influences the structure of BEC) and the Q-dependence 
of BEC represented by correlation function CziQ)- The 
dynamical source of coherence is identified in our case 
with the existence of a constant external term P in the 
Langevin equation. Its influence turns out to be identical 
with the one obtained before in the QS approach [4] and 
is described by eq. (12). Its action is to order phases 
of the produced secondaries. Therefore for P — > oo we 
have all phases aligned in the same way and C2(Q) = 1. 
This is because both here and in [4] the coherence has 
already been introduced on the level of a hadronizing 
source, as property of fields (in [4]) or operators describ- 
ing produced particles. Dividing instead the hadronizing 
source itself into coherent and chaotic subsources leads 
to results obtained in [5] and given by eq. (16). The con- 
troversy between results given by [4] and [5] is therefore 
explained: both approaches are right, one should only 
remember that they use different descriptions of the no- 
tion of coherence. It is therefore up to the experiment to 
decide which proposition is followed by nature: the sim- 
pler formula (16) or rather the more involved (5) together 
with (12). From Fig. 2 one can see that differences be- 
tween both forms are clearly visible, especially for larger 



values of coherence a, i.e., for lower chaoticity parameter 
P- 
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FIG. 1. Shapes of C2(Q) as given by eq. (15) for different 
choices of smearing functions £l(q) (cf. text for details). 




Q [GeV] 

FIG. 2. Shapes of C2(Q) as given by eq. (15) - upper panel 
and for the truncated version of (15) (without the middle 
term) - lower panel. Gaussian shape of Q(q) was used in both 
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From our presentation it is also clear that the form 
of C2 reflects distributions of the space-time separation 
between the two observed particles rather than the dis- 
tribution of their separate production points [18] (i.e., 
it is Fourier transform of two-particle density profile of 
the hadronizing source, p(ri,r 2 ) = p(r\ — r 2 ), without 
approximating it by the product of single-particle densi- 
ties, as in [3]). 

Finally, we would like to stress that our discussion is so 
far limited to only a single type of secondaries being pro- 
duced. It is also aimed at a description of hadronization 
understood as kinetic freeze-out in some more detailed 
approaches. So far we were not interested in the other 
(highly model dependent) details of the particle produc- 
tion process. This is enough to obtain our general goals, 
i.e., to explain the possible dynamical origin of coherence 
in BEC and the origin of the specific shape of the correla- 
tion C 2 {Q) functions as seen from the QFT perspective. 
Actually, our source of coherence should be regarded as 
being only one possibility [3,4], the others were discussed 
in detail in [5]. It is then plausible that in general de- 
scription of the BEC effect they should be somehow com- 
bined, especially if experimental data would indicate such 
necessity. But to do so our approach should first be gen- 
eralized to allowing, as is the case in [5], for production 
of different types of secondaries and allow also for res- 
onance production and final state interactions (both of 
strong and Coulomb origin). This is, however, outside 
the scope of the present paper. 
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